Photon distribution amplitudes and light-cone wave functions in chiral quark models 
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The leading- and higher-twist distribution amplitudes and light-cone wave functions of real and 
virtual photons are analyzed in chiral quark models. The calculations are performed in the nonlocal 
quark model based on the instanton picture of QCD vacuum, as well as in the spectral quark model 
and the Nambu-Jona-Lasinio model with the Pauli-Villars regulator, which both treat interaction 
of quarks with external fields locally. We find that in all considered models the leading-twist distri- 
bution amplitudes of the real photon defined at the quark-model momentum scale are constant or 
remarkably close to the constant in the x variable, thus are far from the asymptotic limit form. The 
QCD evolution to higher momentum scales is necessary and we carry it out at the leading order 
of the perturbative theory for the leading-twist amplitudes. We provide estimates for the magnetic 
susceptibility of the quark condensate Xm and the coupling /a^., which in the nonlocal model turn 
out to be close to the estimates from QCD sum rules. We find the higher-twist distribution ampli- 
tudes at the quark model scale and compare them to the Wandzura-Wilczek estimates. In addition, 
in the spectral model we evaluate the distribution amplitudes and light-cone wave functions of the 
^ p- meson. 
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I , I. INTRODUCTION 

Ph. 
(D 

• Investigations of hard exclusive processes are essential for our understanding of the internal quark-giuon dynamics 
of hadrons. Theoretically, such studies are based on the assumption of factorization of dynamics at long and short 
distances. The short-distance physics is well elaborated by perturbative methods of QCD and depends on particular 
hard subprocesses. The long-distance dynamics is essentially nonperturbative and within the factorization formalism 
becomes parametrized in terms of hadronic distribution amplitudes (DAs) or their transverse- momentum unintegrated 
generalizations, the light-cone wave functions (LCWF)^. These nonperturbative quantities are universal and are 
defined as vacuum-to-hadron matrix elements of particular nonlocal light-cone quark or quark-gluon operators. The 
evolution of DAs at sufhciently large virtuality is controlled by the rcnormalization scale dependence of the quark 
bilinear operators within the QCD perturbation theory. For leading-order DAs this dependence is governed by 
simple QCD evolution equations of the Efremov-Radyushkin-Brodsky-Lepage (ER-BL) type When 
the normalization scale goes to infinity the DAs reach an ultraviolet fixed point and are uniquely determined by 
perturbative QCD. However, the derivation of the DAs themselves at an initial scale /Xq from first principles is a 
nonperturbative problem and remains a serious challenge. Moreover, at experimentally achievable energies it is 
likely that the normalization scale is low and the nonperturbative effects are essential. It is expected that at such 
low scales the shape of DAs can differ essentially from their asymptotic forms and at present can not be strictly 
predicted from first principles. Nevertheless, there are several well-established methods to get information on genuinely 
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nonperturbative quantities such as hadronic distribution amplitudes at low normalization scale: the QCD sum rules 
method fzL IsllgllTo llTil . the relativistic quark models ll2i|. th e instanton Hquid model [13, QIIE EE Il3 1 the effective 
chiral quark models [iH IT^ , or transverse lattice QCD |2(ll| . First and second moments of pion DAs have also been 
evaluated on Euclidean lattices 

A special class of processes involves either virtual or real photons which provide ideal tools for probing the hadronic 
structure in deep-inelastic and hard exclusive scattering experiments. The reason is that the photon is not an 
eigenstate of QCD but a superposition of the Ui gauge boson and quark-gluon configurations which are suppressed 
by the electromagnetic coupling. While a photon is normally considered a structureless particle, it can fluctuate into 
a charged lepton- and quark-pair states, which can be revealed through interactions with a highly virtual photon. 
Thus, the photon DAs have both electromagnetic and hadronic components. The electromagnetic component can 
be calculated within QED, while the hadronic part must be analyzed with non-perturbative methods of strong 
interactions. 

Despite the wealth of theoretical studies not much experimental information is available to constrain various theo- 
retical predictions (for a recent review see |23)- The hadronic DAs and LCWFs can be measured from hard exclusive 
processes at some high momentum scale. For example the pion and photon DAs may be extracted from hard dijet 
production by incident pions ,25i . t26i] and real photons |27l l2q , correspondingly. The measurements of the photon DA 
and LCWF are based on the method of diffractive dissociation applied already to the study of the pion distributions. 
A differential measurement of the pion DA and LCWF was performed by the Fermilab E791 collaboration by studying 
the diffractive dissociation of high momentum pions into two jets 29]. The recent measurements at HERA provided 
the first evidence that diffractive dissociation of particles can be reliably used to measure the photon DAs. Hence, 
the electromagnetic component of the photon LCWF may be extracted from the study of the exclusive ep epfX^ ^~ 
photoproduction process 13(1 and the hadronic component from the exclusive processes of diffractive photo- or electro- 
production of two pions 131] . This may be considered as a special case of the photon dissociation to dijets when each 
jet consists of one pion. 

The aim of the present work is to study the photon DAs and LCWFs of leading and higher twists at a low- 
momentum renormalization scale in a variety of chiral quark models: the gauged non-local chiral quark model based 
on the instanton picture of QCD vacuum [sS. IsA .34] , the spectral quark model , and the Nambu- Jona-Lasinio 
model with the Pauli-Villars regulator. There are a number of advantages of the chiral model approach in comparison 
with other existing approaches. Besides the incorporation of spontaneous chiral symmetry breaking calculations can 
directly be carried out in the Minkowski space. Actually, the corresponding calculations are manifestly covariant, 
so the subtleties of the light-cone quantization do not arise. Within non-local models an effective resummation of 
quark-gluon condensates of growing dimension is performed. 

Although our results are obviously model dependent, we have theoretical control on several properties that should be 
satisfied a priori^ with the symmetry requirements being the most stringent restriction. The resulting Ward-Takahashi 
identities are properly incorporated and their role in quark models is very important when computing the light-cone 
properties. We stress that in our approach we depart from standard parameterizations and do undertake a genuine 
dynamical calculation. Analogous results for the pion have been obtained earlier in the instanton model |l3l UtI l36j ] , 
in the Nambu-Jona-Lasinio model 19] , in the spectral quark model (35i] . and most recently in the large- A^c Regge 
model 113. 

The paper is organized as follows: in Sect. II we give definitions of the photon DAs in different channels corre- 
sponding to the twist expansion up to the twist-4 level. Necessary information about the chiral quark models that 
we are going to use in the paper is given in Sect. HI. In particular, the corresponding normalization constants are 
evaluated. Section IV contains the formalism of quark model calculations of the photon LCWF and DA and a proper 
identification of the low-energy matrix elements is assessed. In Sect. V we give the explicit expressions for the 
leading-twist photon and p-meson distributions, while the perturbative QCD evolution of the resulting leading-twist 
distributions to the relevant physical scales is carried out in Sect. VI. In Section VII the calculation of the higher 
twist components of the photon distributions is presented. Some technical details of calculations, in particular the 
treatment of regularization and the QCD evolution, are given in Appendices fXl- O 



II. DEFINITIONS AND NOTATIONS 

The distribution amplitudes for the virtual photon are defined via the matrix elements of quark- antiquark bilinear 
operators taken between the vacuum and the one-photon state ]7^('7)) of momentum q and polarization vector e^^'' ^ . 



^ Our definitions and notation follow closely the works of Braun, Ball and coauthors l3Sll39l . 
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It is assumed that the quark and antiquark are separated by the distance 2z and the Ught-Uke Umit ^ is taken at 
a fixed scalar product q ■ z. We denote > for space-Uke vectors and < for the time-hke vectors. The photon 
polarization is always perpendicular to q, thus we have e^'*'-' -9 = 0, while for the case of the real photon one has in 
addition the condition e^'''^ -2 = constraining the photon polarization to the two transverse directions. 

Following Refs. [H [H E^, 113 , we use the light-cone expansion of the matrix elements in order to define the 
invariant amplitudes up to the twist-4 accuracy (z^ terms are neglected) 

(0|g(z)a^,[z, ~z]g(-z)|7^(g)) - (0 \qq\ 0) /l^ (<z') | - g^e^^)) Xm dxe'^i'^ At(^x, q^) + 



+ - 



{z ■ q)- 



■ {qf_,z^ - z^qy) 



dxe'^^'-^'BTix^q^) + 

z ■ q 



j\xe'^''-'CT{x,q^)\, 



(1) 



{0\q{zh^[z,^z]q{-z)\j\q)) = ej^^fl^ {q' 



9m 



9m- 



q ■ z 



dxe'^'i-^Bv{x,q^) + z^ 



q - z f 



(g2 



dxe'^'i-'-Av{x,q^)+ 







ig(z)7^75[z,-z]g(-z)|7^(<7)) =ej37/;('?')eM-/34^^9"^'' [' dxe'^''-'V{x,q^), 

Jo 



(2) 
(3) 



where (0 \qq\ 0) is the quark condensate, Xm is the magnetic susceptibility of the quark condensate, and /a^ is related to 
the first moment of the magnetic susceptibility. The symbol [—z,z] in the matrix elements denotes the path-ordered 
gauge link (Wilson line) for the gluon fields between the points — z and z. In the light-cone gauge, A{z) ■ z — 0, 
assumed throughout our calculations, we have [z, — z] = 1 and the Wilson lines may be dropped. The integration 
variable x corresponds to the momentum fraction carried by the quark and ^ = 2a; — 1 for the short-hand notation. 
The electric charge of the quark is denoted by Cq. For a real photon, due to condition e*^^' • z = 0, only four structures 
corresponding to the invariant amplitudes At, Ct, By, and V survive. In the present work we do not consider the 
twist-3 three-particle DAs which involve the gluon fields, integrated out in effective quark models. 
The corresponding decay constants and form factors are defined as 



{0\q{OWMO)\l\q)) = (0 \qq\ 0) Xm/1^ (q^) [e^^U. ~ qA^^ 
(0|g(0)7M'?(0)l7'(9)) = e,/3,/I, {q') el'\ 



dZf3^ 



\q{z)j^t75[z, -z]g(-z)|7^(g)) 



(4) 
(5) 
(6) 



The vector constant at zero virtuality, f]_^ [q^ ~^)^ ^^""^ conservation of the vector current. The 

nonperturbative constants Xm and /3-y provide natural mass scales for the tensor and vector components of the 
photon DAs. In QCD the matrix elements of the relevant local operators are defined as 



(0 \qa^uq\ ^)f = ^qXm (0 \qq\ 0) i^^^, 



qgG 



(7) 
(8) 



where F^j^^, (G^j/) is the field-strength tensor of the external electromagnetic (gluon) field. Da is the covariant derivative, 
and the index F indicates that the vacuum expectation values are taken in the QCD vacuum in the presence of the 
background field . QCD predicts the scale dependence of the quark condensate (0 \qq\ 0), its magnetic susceptibility 
Xmj and /3.y, which at the leading logarithmic approximation evolve as 



(0|gg|0)| =L-^'/^ (0|gg|0)| 



Xr, 



(9) 



where L — as (/i^) /ctg (/ig) is the evolution ratio, b = (llA^c — 2ny) /3, 7^^ = — 3Cf is the anomalous dimension 
of the quark condensate, 70 = Cp is the anomalous dimension of the chiral-odd local operator of leading- twist. 
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7/ = 3Ca — Cf/3, with Cp = 4/3 and Ca = 3 for Nc = 3. By convention, the constants are chosen in such a way 
that the invariant functions are normahzed by conditions 

/ dxAT{x,q'^) = l, [ dxBT{x,q'^)^ [ dxCT{x,q'^) = 0, (10) 
Jo Jo Jo 

/ dxAv{x,q^) = / dxBvix,q^) = 1, / dxCv{x,q^) = 0, 
Jo Jo Jo 



[ dxV{x,q'^) 
Jo 



1. 



In order to define the photon DAs of definite twist it is convenient to introduce the hght-hke vector (p^ = O) , 
such that 



where the Ught-hke vector 



p^ = q^- (11) 



n, = ^ (12) 
p - z 



is normalized by the condition n ■ p = 1. Then the photon polarization vector e^^^ is decomposed into projections 
onto the two light-like vectors and the transverse plane as 

eW = (eW.n)p,+ (e(^).p)n, + ei^). (13) 

Applying the useful relations 

z-q = z-p, e'-^^ ■ p = (e'-^^ ■ nj , (14) 
it is possible to express the above definitions in the form 

{0\q{z)a^.q{-z)\^HQ)) = ie,{qq)fU [q^) | (e^^. - e^^lp^) Xm dxe'^''-^cj>^^{x, q^)+ 

+ {Pt^n, - n^Pu) (e(^)n) dxe^^'^Vf (a;, q^) + (e^L^^. - n^e^^) dxe'««-/iW(x, q^)^ , (15) 

{0\q{z)j^q{-z)\^\q)) = e,h,fl^ {q^) {e^^) . n) ^-^^j\xe'^'i-^ c^\\,{x,q^)+ 
+e^ll dxe'^'i-^i^^l]{x, q^) + (e(^) • n) dxe'^'i-^h'>^\x, q^)^ , (16) 

(0|g(^)7^75«(-^)|7^(9)) = e W^^^P"^'' rfxe*««- V^'') {x, q^). (17) 

^0 

The DAs and (/iiy^ are of twist-2, ■i/)^* are of twist-3, and h^^'^^ are of twist-4'^. In the case of the real photon 
one has e^"^^ • n = and only <j)_L-~j, tp'f}^, V't"^ ^'^d are coupled. The leading twist functions (j)_\_^{x,q^) and 
(/)||^(a:, g^) describe the photon distribution amplitudes that a virtual photon with momentum q dissociates into a 
quark- antiquark pair at small transverse separation. 



^ This dimensional twist-counting refers to the "dynamical" twist of a matrix-element, as opposed to the "geometric" twist of a (local) 
operator. 
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Comparing H15() - (|17|l with the hght-cone expansions (QJ-© one finds the following relations between the invariant 
amplitudes and the DAs 

0-L7 (a;, q^) = At {x, q^) , ^7 ^ {x, q^) = q^X^Ar {x, q^) + Bt {x, q^) + Ct {x, q^) , 

/i^') {x, q^) = ^qhmAr {x, q^) + Ct {x, q^) , 

011^ (x, g2) ^ Av {x, q^) , V^^") {x, g2) = Bv [x, q^) , 

h(;^ {x, q^) = \q^!l2^^Av (x, q^) q^By {x, q^) + Cy (x, q^) , 

(18) 

Within the quark models it is more convenient to calculate first the invariant amplitudes and then by using the above 
relations the DAs themselves. 

III. CHIRAL QUARK MODELS 
A. Instanton-motivated nonlocal chiral quark model 

The analysis in chiral quark models is carried at the one-quark- loop level, i.e. at leading order in the number of 
colors, Nc. We consider the strict chiral limit assuming zero current quark masses, t nn.d = 0. Our basic expressions 
are derived for the case of the non-local chiral quark model, discussed in detail Ref. [ssL 142. 11^ 143. l45l | . The expressions 
for other models analyzed in this paper can be formally obtained from the general formulas of the non-local model. 

The non-local quark models are inspired by the underlying QCD-based models, such as the instanton- vacuum model 
the Schwinger-Dyson resummation of the rainbow diagrams (50j, and some others, which all may be 
cast in the form of non-local quark dynamics. The models develop, via spontaneous chiral symmetry breaking, the 
quark mass depending on the quark virtuality. The quark propagator has the form (following the instanton liquid 
model we set for simplicity the quark wave- function renormalization to unity, Z{p) = 1) 

p — M (p) + le 

where the dynamical quark mass M(p) is ex pres sed via the nonlocal (in coordinate space) function / (p) defining the 
nonlocal properties of the QCD vacuum [HH l52l IH^ . We have to emphasize the essential difference of this expression 
compared to the usually used perturbative expression, where the effects of the spontaneous violation of the chiral 
symmetry are not taken into account and thus M{p) = 0. 
Throughout the paper we use the notations^ 

D{k)^e + M{k), mW {k,k')= ^^^^^l^f\ (20) 



fc± = fc±k fci = (fc+-fc-)-^^i^#^, (21) 

z q 



and for any function F one defines F± — F {k± ) . 

The photon-quark coupling for the incoming quark of momentum k and the outgoing quark of momentum fc' — k + q 
is equal to ,5^ 

r^(fc, k') = 7. - (fc + k'),Mli' + (j, ~ ^) By {q') , (22) 
where the contribution of the rescattering in the p meson channel is taken into account by the factor |54j 
j.^ (^2^ _ Gv^ g,, f d^k fXfX , , , 2,2 , ,.,2(1),. . , 4,2/-/(^Hfc-,fc+) 



M+M_ +k+-k_- -kl [1 + M2(i) (fc+, fc_) j + -kl — 



(23) 



* In following the integrals over the momentum are calculated by transforming the integration variables into the Euclidean space, (fc'' — > ik4. 
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The polarization operator in the vector channel is given by 



J^{q^) = -2N, j ■ 



^4 D,D_ 



M+M^ + k+ ■ - -k 
o 



(24) 



{2ny 

Importantly, the vertex 1)22(1 satisfies the Ward-Takahashi identity for the vector current, 

q^,T^'{k,k')=S-^{k)-S-\k'). (25) 

The straightforward check uses the fact that By (0) = identically. The inclusion of dressed vertices is crucial for 
the consistency of calculations in the non-local quark models. 

As argued in (HJ] (see also [s^.lssjl'). it is convenient to work with the mass function corresponding to the instanton 
field taken in the axial gauge, where Pexponential factor over nonperturbative (instanton) gluon field equals unity. In 
this gauge at large space-like momenta M{p) must decrease faster than any inverse power of p^, e.g., as an exponential. 
Below we assume for simplicity a universal gaussian form of the two nonlocal functions, 

/(p) = /v(p)=exp(^-^), (26) 

with p denoting the Euclidean momentum. Also in this work we do not consider for the nonlocal model the time-like 
region close to the resonances, such as the p-meson. As the model parameters we take 

Mo = 240 MeV, A = 1110 MeV, (27) 

from the work |56j| where they were fixed by the requirement to reproduce the chiral limit of the pion weak decay 
constant /^'^'''^' = 86 MeV [53 and the contribution of the hadronic vacuum polarization to the anomalous magnetic 
moment of the muon. As found in 56], the latter is very sensitive to the value of the mass parameter Mq, preferring 
lower numbers. For the vector coupling we take Gy — —8.7 GeV~^. In Appendix IbI we give further details on the 
properties of the Gaussian shape of nonlocality in momentum and a-representations. 

Let us also mention restriction on the nonlocal models coming from consistency with the operator product expansion 
in QCD. It was first found in IJl that in order to guarantee correct high momenta expansion of hadronic form factors 
and other quantities only vertices depending on quark virtualities k^ and k'"^ may be used. The possible forms are 
f{k'^)f (fc'^), as it is predicted by instanton model, or more general form /(^-^j^)- The vertices containing the linear 
combination of quark momenta, such as fi ^t^ )' l^ad to exponentially growing hadronic form factors at large external 
momenta. The same requirement restricts possible prescriptions for the usage of the flavor P-exponents gauged with 
respect to external currents. We follow the prescription suggested by Mandelstam and first used within the effective 
models in ■42]. On the other hand, the prescription used in 32] leads to a linear combination of quark momenta 
and results in exponentially growing asymptotics of hadronic form factors, inconsistent with the operator product 
expansion. 

Now we are ready to compute the quantities referring to the quark bilinear matrix elements of Eq. (|1I3|) . The 
normalization factors are defined as (we denote u — k'^)^ 

_ Nc f du u {M (u) - uM' (u)) 



V 47r2 D{u) 

The above equations with parameters H27|) yield 

\ inst 



{qq) J 4^2 ^2 ' (29) 



{Omor^' = -(0.214 GeV)^ xrL. =4.32 GeV-2, (31) 
/i'^^t]^^ = -0.0073 GeV^. 



^ The expression for the magnetic susceptibility has been obtained in l54ll58l . 
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In Ref. [SJI the renormalization scale typical for the instanton fluctuations was estimated as /iinst ~ 0.53 GeV. The 
estimate is made by means of comparing the value of the quark condensate with that found in QCD sum rules at 

standard renormahzation point 1 GeV. With Aqcd for three flavors being Aq'^^^'^^ — 312 MeV |63| the 

evolution ratio in from the scale /iinst to /ip equals L = 2.17. Rescaling the values obtained in the non-local quark 
model to the standard renormalization point fiQ yields 



(0|w|0) 



mst 
37 



= -(0.24 GeV)^ X. 
= -0.0035 GeV^. 



inst I 



2.73 GeV" 



(32) 



These numbers are in a rather good agreement with the estimates obtained from the QGD sum rules supplemented 
by the vector meson dominance technique, which provide j39| | 



{0\qq\0) 



QCDsr 



j-QCDsr 
/37 



Mo 



(0.24 ± 0.02) GeV", x 
-(0.0039 ± 0.0020) GeV^ 



QCDsr I _ 



(3.15 ±0.3) GeV" 



Let us also mention that the vector-meson dominance model predicts |3 

2 



,vmd 



/•vmd 



where rrin 



Xm"" = — = 3.37 GeV-^ f,^ -- 
770 MeV, fp = 215 MeV, and >?3p = 0.10 ± 0.05 



(0.0046 ±0.0020) GeV^ 



(33) 



(34) 




FIG. 1: Diagrammatic representation of the photon-to-current transition. 



The calculations of the form factors consist of the evaluation of the one-loop diagram of Fig.^ where one 

vertex is the photon-quark coupling e^F {k, k') ■ e^^^ (q), and the other vertex involves the operators ciafj, 7p, or 7,^75. 
The photon DAs, discussed later on, are obtained from these one-loop diag rams with the light-cone quark momentum 
projection fixed. The derivation of the form factors is straightforward pj] and yields 



/I7 



d'^k 1 



Xm(0|gg|0)7 4^4 



M4 



k ■ q , 



Af_) 



(1 



Bv {e)fxn 



(35) 



r d*k 1 






/ 47r4 Z)+D_ 





+ 3 



sub 



{l + By (g2)/r/r)-^fciA/«(fc+,fc_) 



(36) 



/I7 il') = 

/37 



r d^k 1 






j Att^ D+D^ 





2 

/l-l ' k " — k I 

+ 3 



[l + Bv {q^)jlf_)-tklNf^'\k^,k.) 



INr 



/7(.^) = l-:^/g^i 



3/37 



Ml 



Ml + Ml' (u_ + Ml) 



+ (+ 



-) 



(37) 
(38) 



The subscript sub means that one needs to subtract the perturbative contribution corresponding to the same expression 
with M{k) set to zero. The results are presented in Fig. |5| We note that except for /°, the form factors drop with 
the characteristic scale of about 1 GeV. 

In the case of the real photon the one-loop calculations provide the normalization constants 



/l7(0)-l, fl,{q^)^0{q^), 

/,r,(o).i, /7(o)^i-|:/^ 



(39) 



du M^ (u) 



8 



The vanishing of the transverse vector part at zero is just the manifestation of the Ward-Takahashi identity. The 
axial constant is not zero because the spontaneous breaking of chiral symmetry is taken into account. 
At large Euclidean q one has 

/i^(g^oo) = — 1, /I^ (g ^ oo) = 0(g-4), (40) 

Xm q 

rv ! A f du UA'P (u) 1 , , 




FIG. 2: Form factors associated with the photon distribution amplitude plotted as functions of the space-like momentum 
squared (solid line for the tensor form factor /l-y, dash-dot line for the vector longitudinal form factor f "^^, dashed line for the 
vector transverse form factor fj^j, and short-dashed line for the axial form factor /°). 



B. Spectral quark model 

The spectral quark model (SQM) is based on the spectral representation of the quark propagator [s^. The 
expressions for one-loop observables in SQM are obtained from the expressions of the preceding sections by replacing 
the quark mass with the spectral mass, M — > w, and then integrating over omega with the spectral density p(w). We 
have 

A= diLjp{uj)A{uj), 
Jc 

where C is a suitably chosen contour in the complex uj plane. In the meson-dominance version of SQM |35| we have 

p{uj) ^ pv{oj) + psi^^), (41) 
11 1 

pv(i^) 



27riw(l-4w2/Af2)5/ 



( ^ 1 48^^(0 ItoIO) 

^^(")--27rz7V,A/|(l-4c.VM|)^/^' ^ ^ 

where My is the mass of the p meson, and ms provides the scale for the scalar spectral density. We note that the quark 
condensate is a model parameter. Matching the SQM predictions to the large- TVc results of the chiral perturbation 
theory in the single resonance exchange approximation allows for the identification [6ll | 



Ms = Mv = mp. 
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w-Complex Plane 



-M/2 



M/2 



FIG. 3: The contour C for the ti;-integration in the spectral quark model. M denotes the p-meson mass. 

where mp is the mass of the p meson. The contour C for the integration over the spectral mass uj is displayed in 
Fig.E 

SQM with the meson dominance (|42|1 generates, by construction, the monopole form of the pion electromagnetic 
form factor [ssj . Interestingly, the model has the feature of the analytic quark confinement |62| . i.e. the quark 
propagator has no poles, only cuts, in the complex momentum plane. Moreo ver, the evaluation of low-energy matrix 
elements in SQM is very simple and leads to numerous results reported in Ref. '35], in particular for the pion light-cone 
wave function and the pion structure functions. 

The results for the normalization factors in SQM are 



SQM I 



^ 2 =3.37GeV-^ /f^^^ 
mi ' 



Nr m 



MSQM 



47r2 6 



= -0.0075 GeV . 



In Ref. [33 it has been argued that the quark model scale corresponding to SQM is very low, /^sqm = 313 MeV. That 
yields a large evolution ratio L — 4.6 between /^sqm and fiQ —1 GeV, and, consequently, we find after evolution 



C^Q^I, =1.37GeV-^ /: 



SQM 

37 



= -0.0018 GeV^ 



much lower values than in the nonlocal model and QCD sum rule. 




0.5 1 1.5 



q [GeV] 



FIG. 4: Same as Fig. |5]for the spectral quark model. 
For the tensor form factor we find 

/I7 ii'^) = ^2 r„2' 
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in full accordance to the vector dominance model. This makes a strong case for the assumed scalar part of the 
meson- dominance of spectral quark model, with the power dg = 5/2 (see Ref. H^)- The tensor form factor satisfies 
the conditions 

Amy 

For the axial form factor we find in local models /"(g^) = 1. Expressions for the vector form factors, /^|| and f^j_ are, 
due to the vacuum subtraction, not analytic and the integration over the x variable has to be carried out numerically. 
We present these form factors for SQM in Fig. 01 noting the same qualitative behavior as for the non-local model in 
Fig.H 



C. Nambu— Jona-Lasinio model 



The expressions for one-loop observables in the Nambu- Jona-Lasinio model (NJL) with the Pauli-Villars regulator 
follow from the expressions for the non-local model. Formally, one replaces in the denominator the square of the 
(constant) mass, -I- A^. Then, in the simplest twice subtracted case one has for any observable A the 

prescription 

A = AiA^ = 0) - A(A2) + A2 "^^^^^^ 



dA2 



The parameters used are M = 280 MeV and A = 871 MeV, which were adjusted to yield F^r — 93 MeV. The used 
NJL model has no explicit coupling to vector mesons (Gy — 0). For the normalization factors we get 



(0|gg|0)i 



NcM 
47r2 



A^ - log 1 



A^\ 



= -(0.230 GeV)^ 



NcM 



47r2 (0 \qq\ 0) 



A2 



^37 I 



NcM^ 
47r2 



A2 



AP + A 



AP + A2 
-logfl 



- log 1 



A^ 

AP 



2.55 GeV~^ 



A^ 

AP 



= -0.0087 GeV^ 



The QCD evolution brings these numbers down, similarly to the case of SQM. 
For the tensor form factor we find 



(A/2 -I- A2) 2(2A2 + 5'o)tanh"^^ - 2V3b5'tanh"^;|^ + qy/Slog (l 



IP 



qVS [(M2 A2)log(l 



A2 



A2] 



S = 4Af 2 + 4A 



2 ' q^, 5o =4M2 + 4A2, 



(44) 



(45) 
(46) 



which satisfies H43|) . Other form factors involve numerical integration. All form factors obtained in the NJL model 
are qualitatively very similar to the case of SQM shown in Fig. 01 



IV. PHOTON LIGHT-CONE WAVE FUNCTIONS IN QUARK MODELS 

We use the light-cone coordinates with the convention a± — oq ± 03. The null vector = (nP , , ri^ , rv^) = 
(l/p+, 0, 0, — l/p+) satisfies the conditions — and n ■ p ^ 1, n ■ a = /p^. Wc also denote 

d'^k dk+dk^d^k±^ _ ~ 
(2^0^ " 2(2^)4 " 

The definition (|15|l is equivalent to (from now on we only keep the leading twist contribution) 

»e,Xm(w)/l^ {q^) {qpe\^^ - qJp^) Mx) = T '^e-^^^^^-^HQ\q{Tn)ac,pq{-Tn)\^\q)) . (47) 
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In the one loop approximation the quark model evaluation of the matrix element of the quark bilinear yields 

Xm(w)/1^ (g') [qpei^^ - fce^-")) - -iN, J dk6{k ■ n - x)Tr[fTa^5(fc)r • e'^^^Sik - q)], (48) 

where the trace is over the Dirac space. The interpretation of this result is clear: the photon DA times a Lorentz 
structure is equal to the quark-loop integral, where the +- component of the momentum of one of the quarks is 
constrained to k+ = xp+ (see Fig. If one does not carry out the transverse momentum integral in Eq. (|48|l . then 
one obtains the photon light- cone wave Junction, $j^-y(a;, fcj^), 

Xm(g<?)/1^ (g') (q0e(^^-qJp^)'^±,{x,k^) = ^iN,j '^^^5{k ■ n - x)Tr[a^pS{k)T ■ e^>^'> S{k - q)]. (49) 

Obviously, we have 

j (i^fcj_$j_-y(x, A:j_) = (/ij_^(x). 

For the twist-3 vector component of the photon light-cone wave function and DA we find from the definition llt)|l . 
with analogous steps as above, the expression 

fc^) + q,n ■ e(^)*2(2;, ki_) ^ ^ J ^^^^n ■ k - x)Tr[j,S{k)T ■ e^Sik - g)], (50) 

while the DAs are ^5^2:) = / d^/cj.^'J^ (a;, /cj.), with i = (||,±). 

The case of the axial twist-3 component is somewhat more complicated, since the definition (|17|1 involves a power 
of the coordinate z in the tensor structure. As a result we find 

ej3,f^ (q^) e^'"'^ei^^q^np^i;\x) = . . (^--i)(0|g(rn)7^75g(-Tn)|7^(g)). (51) 

The one-quark loop evaluation yields 

4iV, 
hi 

We may introduce the derivative of the distribution, d / dx ipif°'\x) , for which we have 

e,f3if^ e^"""'ei^kc.np-^^i;\x) = ^J dk5{n ■ k ~ x)T^[r-i,S{k)T ■ e^^^ S{k - q)]. (53) 

The evaluation of this quantity is similar to the evaluation of (j)^ or -i/;'^"-' . The axial light-cone wave function is obtained 
from Eq. if the fcj, integration is left out, in analogy to Eqs. H49I5U|I . 

V. PHOTON LIGHT-CONE WAVE FUNCTIONS IN CHIRAL QUARK MODELS 

A. Leading-twist photon DAs 

The distribution amplitudes of the real photon calculated in the instanton model in the chiral limit may be cast in 
a closed form. It is convenient to introduce notations for the integration variables^ 

u+ = u — iXx, U- = u + iXx, M± = M {u±) , D± = D {u±) , x = \ — x. (54) 

By using the definitions of Sect. IV one gets the twist-2 DA for the cr^i, structure (Fig. O 



.h^r; {q') e^''''^ei^^q^np^{;^\x) = / d~ke{n ■ k - x)Tiirj,S{k)T ■ e^^'>S{k - q)]. (52) 



(t>A^^{x,q^ = Q) = 1 e(xx) I du-^^^-j du I — n n 

Xm 47r 







D (u) Jo 7_oo 2^ D+D 



(55) 



' Some details of calculation method first developed in llfHSRl and the relevant expressions for the integrals are given in Appendix 1X1 and 
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d> (x) 



0,0 



0,0 0,2 0,4 0,6 0,8 1,0 



FIG. 5: The leading-twist photon distribution amplitude in the tensor channel, (^^^(x,q^ — 0), evaluated at the quark model 
scale in the non-local quark model - solid line, in the spectral quark model and the NJL model {= 1) - dot dashed line, its 
asymptotic form - dashed line, and the result of the local approximation to the instanton model - dotted line. 



Within the non-local model we can find from H55|) the following light-cone wave function: 



Xm 47r2 



e(xa;)- 



M (k1 



2^ D+D- 



M(i) {u+,U-) 



(56) 



It exhibits an exponential decay law at large k^ , which is a consequence of the nonlocality of the model |54| and is in 
accordance with the conclusion made in that finiteness of all kj_ moments of the quark distributions guarantees 
the exponential fall-off of the cross sections. 

In the vector channel the leading twist structure for the real photon < 
zero itself and has form 

hi (^'9^ = 0) = Q{xx) 



(x, — O) is decoupled but the DA is not 



(57) 



independently on the nonlocality shape. 

The results for SQM and NJL models are obtained from (|55|l by taking the local limit (all derivatives are zero 
M' = 0, M(i) = 0, etc. ). This gives the result 



,g2 = 0) =e(xa;). 



(58) 



for i =_L, II . These leading twist DAs calculated within the instanton and local models are shown in Fig. [Sjalong with 
its asymptotic form defined at very large scale 



/asympt / \ — 



In Fig. |S1 there is also given the result obtained in |l6l | in gauge noninvariant approximation 

d\{l-x)M++xM_ 



A; 



approx 



du 



2tt 



D+D^ 



(59) 



(60) 



47r2 



du 



u [M (it) - \uM' {u)) 
DHu) 



when the full vector vertex (|22|l is substituted by the local part. This approximation resulted in a violation of the 
proper normalization as implied by the Ward-Takahashi identities. 

The dependence of the leading twist photon DA {x,q^^ on the photon virtual! ty in space- like and time- like 
regions is demonstrated in Fig. |B1 Both leading twist DAs are normalized by unity 



07 = 1- 



(61) 
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FIG. 6: Dependence of the twist-2 tensor component of the 
photon DA on transverse momentum squared {cf = 0.25 
GeV^ sohd line, = GeV^ dashed hue, = -0.09 GeV^ 
short-dashed hue, asymptotic DA - dotted line). Nonlocal 
model at the quark model scale. 



1,50 




0,75 



0,50 - 
0,25 - 

0,00 I 1 1 1 1 1 

0,0 0,2 0,4 0,6 0,8 1,0 

X 



FIG. 7: Same as Fig. |^for the twist-2 vector component of 
the photon DA. 



At very high space-like photon virtuality the photon DA approaches the (5-function distribution concentrated at the 
edge points of the x interval. It corresponds to the configuaration when all momenta pass through quark or anti- 
quark line, with other line being soft. At time-like virtualities approaching the p- meson pole the photon DA will be 
concentrated in the middle of x interval. 



B. Spectral quark model 

The expressions for SQM are obtained from the formulas of the preceding sections with the formal replacement 
^{p) ~^ and supplying the loop integrals with the addition spectral integration p{uj)duj, as explained in 
Sect. III.B. As the result, we find the following expression for the leading-twist component of the light-cone wave 
function for the real photon in the tensor channel: 

"^^^^"'^^^^ ^^(l + JiMg)^/^ ' ^''^ 

Note the power-law fall-off at large transverse momenta, ^^{x,k±) ^ 1/k^. In cross section this leads to tails of the 
form ~ 1/kY'- 

For the virtual photon of virtuality q^, the tensor component of the light-cone wave function has the form 

6(1 + ^. 



i>±r = : \. ..5/2 ' (63) 



mlll + 4 '- 



We also find 

2TOp(m2 + q^){3ml + q"^ + 12kl) 

{m?p + 4fe2^)3/2(TO2 -|_ g2 _|_ 4fc2^)2 



"'0 



(64) 



The transverse-momentum rms is 



TO 



2 



0j_^-(x) 2 

For the real photon one has the estimate {k\) ~ (544 MeV)^. For the instanton model one gets 
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For the real photon DA one finds the constant form (|58|l . while for the virtual case we find the following normalized 
DA: 



1 + ^ 



I I | g^a;(l-z) 



3/2 ■ 



The corresponding expressions for the p meson are obtained by taking the limit ^'^pi which for the DA yields 

Mx) = S{x - i) (67) 



This highly singular behavior is washed out by the QCD evolution, see Sect. IV Dl 

The leading twist photon amplitude in the vector channel, (f)\\^{x), is of the form given in Eq. (|58|l . i.e. constant. 



C. NJL model 



In the NJL model one obtains the same leading-twist DAs for the real photon as Eq. (|58|l . i.e. {x,q^ = O) = 
Q{xx). Extension to the p- meson pole, however, is problematic due to appearance of the quark production threshold. 
For this reason we do not perform this analysis. Note that SQM is free of this problem, as no quark poles are present 
in this model - a feature sometimes referred to as analytic confinement. 



D. QCD evolution 



The results of the previous sections referred to the low-energy scale of the quark model, fi. In order to relate to 
results at higher scales, the QCD evolution is necessary. Within the effective model approach the method has been 
described in detail in Ref. [13 and is similar to the case of parton distribution evolution [sM l6^ . For the twist-2 
photon and p-meson DAs the leading-order QCD evolution is made on the basis of the Gegenbauer polynomials. The 
basic logic here is that the quark models provide the initial condition for the QCD evolution, now written as /i^). 
The evolved leading-twist distribution amplitudes read 65] 

^\x,q^) = 0L(2;)^'r=oC,t(2x- l)aUg'), (68) 

with 

(t'Isix) = 4>Isix) = 0a7(a;) = Qxx, 

denoting the Gegenbauer polynomials, the prime indicating that only even n enter the sum, and A = 3/2 for T, 
V, and AV. The Gegenbauer coefficients ajj evolve with the scale. The corresponding formulas for the leading-order 
QCD evolution are given in Add. lUl We note that the exponent in Eq. (|C1|) contains the difference (7,^^ — 7q), such 
as to make the zeroth moment (the norm) of the DAs constant. 

Our methods of determining the a priori unknown quark-model scale ^ have been presented in for the case of 
the nonlocal model and in Ref. for the local models. 

The numerical analysis of the evolution of the tensor DA for the photon and the p-meson proceeds exactly as written 
in the above formulas, with the series approximated by finite sums over ri. Numerically, about 100 terms are needed 
to achieve the accuracy in the presented figures. The results for the twist-2 photon DAs are shown in Figs.|S|and|^ 

The results for SQM are shown in Fig.^| We show the evolution of the DAs for the real photon, for the p meson, 
as well as for virtual photon. We note that the evolution radically changes the character of the curves. For the real 
photon one starts with the fiat DA of Eq. (|58|l . exactly as for the case of the pion in the local model 1^, with the 
evolution providing vanishing of the amplitude at the end-points and gradually approaching the asymptotic limit. 
So it means that all moments of the DAs evolve as continuous function, but the DA itself at some points (edge 
points in the discussed case) evolve discontinuously. For the p-mcson case the result is even more dramatic, with the 
initial condition having the singular form of Eq. H67f) . evolving to broader regular functions and also tending at large 
evolution scales Q to the asymptotic limit. We note that this evolution is very slow. The lower part of Fig. I lUI shows 
the DA for the virtual photon with two selected virtualities: time-like q'^ — — m^/2 and space-like = Sm^. We have 
to note that because of low initial evolution scale the QCD evolution carried out at two-loop level would give more 
realistic results. 

The results for the real photon in the NJL model are the same as in SQM. 
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VI. HIGHER TWIST COMPONENTS 
A. Non-leading twist photon DAs in the non-local chiral quark model 

For the real photon the DAs of higher twists are following: the twist-3 transversal DA (7^ structure) (Fig. Ill|l 



5{x) +5{x)^ 26 [xx) + /g-^ ^ / du 



00 /'OO 



d\ (M_ - M+Y 
2^ Di D_ 



the twist-3 DA (7^75 structure) (Fig. EJ) 



1 Nr f°° 

du 



j{a) 4^2 



/37 



(0)^ 



dA xMt + xMl 
2^ Di 



(69) 
(70) 

(71) 
(72) 



r=l, 0.5, 0.25, 0.1, 0.05 



r=l, 0.5, 0.25, 0.1, 0.05 




0.2 0.4 




r=l, 0.5, 0.25, 0.1, 0.05 




0.2 0.4 0.6 0.8 1 



FIG. 8: The LO ERBL evolution of the nonlocal model predictions for the leading-twist tensor photon DA (j)jj^{x,q'^). Top 
left: the real photon DA, top right: the virtual photon at = 0.25 GeV^; bottom: the virtual photon at = —0.09 GeV^. 
The dashed lines show the asymptotic DA, 6x{l — x). Initial conditions, indicated by dotted lines, are evaluated in the nonlocal 
quark model at the initial scale /i'"'^* = 530 MeV. The solid lines correspond to evolved DA'a at scales Q — 1, 2.4, 10, and 
1000 GeV. The corresponding values of the evolution ratio r are given in the figures. The appearance of tiny wiggles is a 
numerical artefact. 
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the twist-4 DAs {cr^^ structure) 



h!^ (x,g2 = 0) = i 

/ dxh^'^ (x, q^) = 0, 
Jo 



5{x) + 5{x) 



JMl^Jo "loo 2^1^^^ 



(Ml - Ml) (Af_ -M+y 



D+D^ 



(73) 



V^W(x,g2=0)=/i(*)(x,9' = 0). (74) 
In the case of Gaussian form factor H26I chosen in the present paper ) the integral in (|73|l may be partiaUy done as 



°° , f°° dX fM+ M\ M„A^ 



A2 



--2nj Q (- ~2n] + {x ^x) 



(75) 



We may compare our results with the Wandzura-Wilczek (WW) approximation which attempts to reconstruct 
non-leading twist result from the leading contributions. The corresponding relations for the photon deduced from the 
p-meson case 0| are 



^(a)WW (^)^2 



du 



(76) 
(77) 
(78) 



From the fact that both the nonlocal model as well as the spectral model predict for the leading twist photon DAs 
an almost constant function, we may predict the nonleading twist DAs as they follow from WW approximation: 



^{a)ww ^ ~2{xlnx + xlnx), 



1-x 



(79) 
(80) 
(81) 



r=l, 0.5, 0.25, 0.1, 0.05 



r=l, 0.5, 0.25, 0.1, 0.05 





0.2 0.4 0.6 0.8 



FIG. 9: The LO ERBL evolution of the nonlocal model predictions for the leading-twist vector DAs of the photon (f)^^_l {x,q'^). 
Left: the real photon DA; right: the virtual photon at = 0.25 GeV^. The dashed lines show the asymptotic DA, 6x{l — x). 
Initial conditions, indicated by dotted lines, are evaluated in the nonlocal quark model at the initial scale /i'"'^' — 530 MeV. 
The solid lines correspond to evolved DA'a at scales Q = 1, 2.4, 10, and 1000 GeV. 
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which should be compared with the exact result obtained in the local models 

^WSQM ^^{S{x)+S{l-x))-e (xx) . (82) 

Taking literally these results one gets that the normalization conditions are violated in the WW approximation. 
Indeed, all ip^^ (x) are normalized to unity. One can fix this problem by shifting the corresponding distributions by 

constants, forcing the norms for ip^"^^^ and ^/;^*'^^ (x) to be zero. The corresponding distributions are shown in 
Figs. El and 1121 

There were attempts to find photon distribution function within next-to-leading conformal expansion with non- 
perturbative coefficients determined from the QCD sum rules |39| |. To the next-to-leading order in the conformal 
expansion the twist-3 photon DA were found as 

(^) = 5 {3e - 1) + tI (15< - 5c.:^) (3 - 30e + 35^) , (83) 



and the twist-4 photon DA 



<^ (-) = ^ (1 - e) i^e - 1) (i + ^< - ^-^^ ) , (84) 



(x) ^ -10 (1 + 2x+) C^^ (^) , (85) 




FIG. 10: Evolution of the real photon DA (top left), p-meson DA (top right), virtual photon at — —'m'p/2 (bottom left), and 
virtual photon at = 2m^ (bottom right). Initial conditions, indicated by dotted lines, are evaluated in the spectral quark 
model at the initial scale p^*^^ = 313 MeV. For the p-meson case the initial condition is given in Eq. I|67|l . The solid lines 
correspond to evolved DA'a at scales Q = 1, 2.4, 10, and 1000 GeV. With the larger the scale the evolved DA is closer to the 
asymptotic form 6x{l — x), plotted with the dashed line. The corresponding values of the evolution ratio r are given in the 
figures. The appearance of tiny wiggles is a numerical artefact. 
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a, ^1 




FIG. 11: Dependence of the regular part (without (5-functions) 
of the twist-3 vector component of the photon DA on transverse 
momentum squared (g^ — 0.25 GeV^ soUd Une, = GeV^ 
dashed Une). The short dashed hue corresponds to the WW 
approximation 11801 and the dotted hue is for the conformal ap- 
proximation 18311 



FIG. 12: Dependence of the twist-3 axial component of the pho- 
ton DA on transverse momentum squared {q^ — 0.25 GeV^ solid 
line, q^ — GeV^ dashed line) 



where the parameters u)}^ , u)^ and >c'^ correspond to matrix-elements of local operators of dimension 6. In the 
vector-dominance approximation they are 



w^ = 3.8±1.8, tj;f = -2.1± 1.0, x+ = 0. 
Note that at low scales the existing QCD sum rule results do not provide reliable results [s^ . 



(86) 



v.(x) 
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FIG. 13: Dependence of the regular part of the 
twist-3 tensor component of the photon DA on 
transverse momentum squared {q^ = 0.25 GeV^ 
solid line, = GeV^ dashed line, q^ = -0.09 
GeV^ short-dashed line). The dotted line is for 
the WW approximation II81II . 



FIG. 14: Dependence of the regular part of the 
twist-4 tensor component of the photon DA on 
transverse momentum squared {q^ — 0.25 GeV^ 
solid line, = GeV^ dashed line) 



The QCD evolution of the higher-twist components of the DAs is a very complicated problem. The reason is that 
unlike the leading-twist case, where the evolution involves the single distribution and leads to simple ERBL equations, 
the evolution of higher twist components couples the quark bilinears to the three-particle quark-gluon components 
p. ISSL l66l| . At the moment there is no manageable procedure we could promptly use, hence we cannot repeat the 
evolution calculation for the quark model prediction for non-leading twist DAs. This important problem is interesting 
in its own right, but it extends outside of the scope of this paper. 
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VII. CONCLUSIONS 

In the present work we have carried out an analysis of the hadronic part of the photon DAs of leading and higher 
twists up to fourth order at low normalization scale within several effective chiral quark models. We have studied the 
DAs of real and virtual photons for both transverse and longitudinal polarizations. Our analysis has been performed 
within two classes of models: the nonlocal model, based on the instanton picture of the QCD vacuum, and local 
models (the spectral quark model and the Nambu-Jona-Lasinio model). Despite the different nature and details of 
the calculations, our analysis has shown a remarkable similarity of the results in both classes of models. By working 
with the photon vertex satisfying the Ward-Takahashi identity we have revised and improved the results for the 
leading- twist photon DA first performed in 16] . We have shown that the calculations based on the nonlocal model 
with a correct implementation of gauge invariance resulted in much flatter results, with the photon DA remaining 
non-zero at the end-points x — 0,1, in qualitative agreement with local quark model calculations. 

The main findings may be summarized as follows. 

(1) The leading- twist DAs of the real photon at low momenta scale are constant in the local models and almost 
constant for the nonlocal model. The independence of the photon DAs on x is expected for the local models. 
They also predict a constant behaviour for the leading-twist pion DA p^. At the same time the nonlocal model 
predicts nontrivial x dependence in the case of the pion |l3l 36] , and shows a constant behaviour for the photon 
independently on the shape of the nonlocal form factor. This is a consequence of the fact that the photon has no 
hadronic form factor in contrary to hadrons which are a genuine bound state of dynamical quarks. 

(2) Furthermore we show that applying the QCD evolution to leading order of the perturbation theory, the leading- 
twist photon DAs become immediately zero at the edge points of the x interval (as in the pion case of Ref. T^) and 
are always wider (flatter) than the asymptotic DA. At the same time the evolution of the moments of DAs is evidently 
continous. This is a typical situation for functional series which are non-uniformly convergent. The role of the QCD 
evolution is very important, as the quark-model results are far from the asymptotic form. The quark model scale 
estimated in Ref. :19] for the local models is very low, about 320 MeV, therefore the evolution is very fast. For the 
instanton model it is estimated as 530 MeV |54ll67j . In short, our predictions may be compared to data only after a 
suitable QCD evolution has been carried out. Further, for more realistic predictions the perturbative evolution with 
a kernel calculated at the two-loop level is certainly desirable. 

(3) Some of the higher-twist photon DAs have delta-function behaviour at the edge points of the x interval. It might 
be interpreted as originating from differentiating the constant leading twist DAs and corresponds to the quark-anti- 
quark configurations of the photon DA when all momentum is carried by one quark (or anti-quark) . This singular 
initial behaviour is washed out during perturbative QCD evolution to higher momentum scales. We also predict the 
behaviour of the leading- and higher-twist DAs for nonzero virtuality of the photon. 

(4) We have drawn our attention to the fact that the equations of motion of perturbative QCD and of the effective 
approaches are different. In the latter case spontaneous breaking of chiral symmetry is taken into account in terms of a 
momentum dependent dynamical quark mass. So our results do not coincide with those obtained in perturbative QCD 
using the equations of motion. In particular, the relations between normalization constants, and DAs for different 
twist (the Wandzura-Wilczek relations) are not satisfied in the effective approach. We question the validity of such 
a method, since it implicitly assumes that current quarks are on the mass shell and free, while our interest is in low 
energy matrix elements. We show that the norm of the tensor amplitude is related to the magnetic susceptibility of 
the quark condensate, and the zero norm in the vector channel is constrained by the vector-current conservation, while 
the axial-vector norm remains unconstrained. The latter is determined by nonperturbative dynamics, and therefore 
its value is model-dependent. 

(5) Finally we compare our results with those obtained from the conformal twist expansion [39l ] using QCD sum rules 
and vector-meson-dominance methods. It seems that the results obtained in the leading and next-to-leading orders 
of the conformal expansion do not converge. We have also considered the Wandzura-Wilczek method of estimation 
of the higher twist DAs assuming a constant distribution in the leading order. 

Acknowledgments 

AD thanks for partial support from the Russian Foundation for Basic Research projects No. 04-02-16445, Scient. 
School grant 4476.2006.2 and the JINR Bogoliubov-Infeld program. This research is supported by the Polish Ministry 
of Education and Science, grants 2P03B 02828 and 2P03B 05925, by the Spanish Ministerio de Asuntos Exteriores and 
the Pohsh Ministry of Education and Science, project 4990/R04/05, by the Spanish DGI and FEDER funds with grant 
no. FIS2005-00810, Junta de Andalucia grant No. FQM-225, and EU RTN Contract CT2002-0311 (EURIDICE). 



20 



APPENDIX A: USEFUL INTEGRALS 



The calculation in the nonlocal models is most involved, hence we show it first. Formulas for the local models 
may then be obtained from the nonlocal model expressions as limiting cases, with SQM involving an extra spectral 
integration. We use the propagator (|19|l and the vertex H22|l . and evaluate the Dirac traces in H49|l for the tensor 
component. The electromagnetic vertex is split into the local part and non-local part, according to Eq. H22|) . We use 
the identity 

2k ■ q = D+ - + Ml - Ml. 

For the local part we find 

j d~kS{k ■ n - i0Tr[a„/35+7 • e^^'^.] = [q^ei^^ - q^ef^jlg + [npei^^ - n^eJ.-'V", (Al) 



with 



Iq = Ai / dk5{k -n ^)- ^ 



D+D. 



(A2) 



2i / dk5{k -n - -S) 



{M+ - MJ) 



1 

1^4 



1 



{M+ - M_)(M2 - Ml +q^{x-x)) 
D+D^ 

The piece involving the nonlocal of the quark-photon interaction vertex has the form 



1 



dkS{k ■ n - ^^)TT[aa,pS{+)S{-)]M[^^_{2k - q) ■ e'-^\ 



(A3) 



All quantities appearing in the above integrals are functions of squares of momenta, k^ — kj^k^ — u and (fc — g)^ ~ 
{k+~q^){k- —q-) — u, with u — k\. Since the transverse and (-1 — ) spaces are factorized, for the clarity of notation we 
introduce K — (fco, k^) and Q — {Qq, Qs)- In order to proceed further it is convenient to pass to the Laplace-transform 
space 



F+ = F(fc2) = 



dae 



G^^G{{k-qf)= d/3e~^(^-Q)'G„(M). 



(A4) 



The terms in Eq. (|A2I) and in other formulas for other observables discussed in this paper have the form of integrals 
over J 5{k ■ n — of functions of fc^, functions of (fc — g)^, or the product of the two. 
Using the definitions ljA4p we find 



d^K5{k ■ n - -0F+ = / d'^K 



d\ r°° 

-oo 27r io 



daFa{u)e 



-aK^ — iXik-n—x) 



where on the way we have carried the gaussian integration over d^K . Similarly. 

j d^KS{k ■ n - \0G- = TT y dzGu+z6{x). 
For integrals involving the product of F-^-G^ we have 



(A5) 



d'K5{k-n - ^S,)F+G- = j d'K j ^ f da f dpFa{u)Gf3{ 



-aK^-l3(K-QY-i\(k-n-x) ^ 



d\ f dadp 



■Fa(u)Gp{u)e 



2t: J{a + 13)- 

Fa(u)Gi3[u)e ^ ^ 



2tt J{a + py 



(A6) 



21 



where A' = A/(a + (3). Since / ^ exp {—iX' {xf3 — xa)) = 6{x(3 — xa), and a > 0, /3 > 0, we verify that the integral 
(|A6I) has the correct support^ Q{xx). 

For the case of real photons, = 0, we have 



(fK6{k ■ n - -^)F+G- / dz 



dX 
2^ 



dad(3Fa{u)G i3(u)e 



— (q+/3)2 — iA(2;/3— xa) 



t: dz 



dX 
2^' 



For the general case of virtual photons one cannot unfold the Laplace transforms. A formula convenient for numerical 
calculations has the form 



d^KS{k-n - -^)F+G- =TT 



(icv ~ ~ _ 2 

a 



(A7) 



For the case on the non-local vertex we encounter integrands which are no longer separable in the and {K — Q)^ 
variables, as they M^^^ of Eq. (|22|l . This quantity may be written as 



m|^L = [ dtM'{tK^ + {l^t){K -Qf +u), 

JO 



where the prime defines the derivative with respect to the argument. Hence the Laplace transform of in the 

variable 7 is —7^(7). 
We introduce 

A = a + /3 + 7, a' = a + 7<, /3' = /3 + 7(1 - i). 
The following integrals are needed in our analysis: 



Ii= j d^K5{k ■ n - i$)F+G_Af|^.Lfc • e^'^^ 



^2 — 



^^^^^F„HG,(u)[-7M;N]e-^^^-^('^'--"'-) 
27r A 

d'^KS{k ■ n ~ x)F+G-M^y [qt^k" - q^k^'] k ■ e^""^ = 



-AXn ■ e('*) 



(A8) 



27r 



Carrying the A integration leads to 5{f3'x — a'x), which together with the positivity of a' and (3' leads to the correct 
support IjASp . 

For the case of real photons we find 



/i -0, 



^2 



(A9) 



dz z 



dX 
2^' 



^ ]\/r\. 



(1) 



APPENDIX B: GAUSSIAN NONLOCALITY 



In this appendix we collect some useful relations valid for Gaussian shape of the nonlocal function f{p) which 
provides the expression for the dynamical mass as (in chiral limit) 

M (/) = Mo exp {~2p^/A^) . (Bl) 

It is important to stress that exponentially decreasing function at large Euclidean momenta is a direct consequence 
of requirements of the gauge invariance and continuity of nonlocal matrix elements ^15- .54j . First of all consider the 
propagator function 

^^^^^ + A/2 exp (-4p2 /A2) ^M^p + exp {-Ap^y^) ' ^^^^ 
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where the second equation is rewritten in the dimensionless variable p = p/Mq, with y = Mq/A. Depending on the 
values of parameters Mq and A the function P may or may not have a real pole. The critical value is obtained from 
equating denominator and its derivative in p^ to zero. Then one gets the region of parameters where the propagator 
has no poles at real values of p"^ as [sj 

y>yo = 1/(2V^) « 0.303. (B3) 

In the other region, y < y^, the real valued pole appears at 

pI = -4KAf2, (B4) 

where the coefficient k monotonically changes between k — 1 obtained in the local limit y — > and n = yfe reached 
at the critical point y yo. 

In the present paper in the nonlocal model we use the parameters (|27|l corresponding to y* = 0.195, where the pole 
(|B4|I appears with k = 1.13 at p§ — —0.29 GeV^. This is a reason why our calculations are given in the region far 
away from the pole. The chosen value of the vector coupling Gy corresponds to the p-meson mass at the somewhat 
low pole position of Mp = 541 MeV. 

The propagator function has also an infinite number of complex valued poles (a similar situation is also encountered 
in the NJL model with the proper-time regularization [b^). This fact makes it difficult to carry out calculations of 
integrals directly in the momentum space. A big advantage is to transform all integrals to the a-representation where 
some analytical calculation may be done [T5I . W(\ . To this end we expand (|B2|) as 



P(.)^^ V I -^^o^expH./A^)y ^^^^ 




and apply the Laplace transfers 



1 1 /I „ (a — A)'^ ^ „ , 



to each term of the expansion. Then one gets the a-representation of the propagator function 



P is) P (.) . 1 + £ [-Mo-(^^;fe/A-)]e (« - 4./A^) . (B6) 

fc=i 



This transformation may be done for all functions. For example, one has transformation 
exp(-./A^) ~ ^ [.Mi{a-{Ak-3)/A^)]'~\ , 

The functions P (a) and P/ (a) are presented in Fig. ^1 Let us note some general properties of the transformed 
functions. At small a one has P (0) = 1 because large s asymptotic P (s ^ 00) = 1/s. In fact for Gaussian function 
f{p^) from ljB6|) and ljB7|l one has P = 1 in the interval a € [0, 4/A^] . At the same time P/ (0) = 0, since at large s the 
function Pf (s) decreases faster than 1/s. More precisely one has for Gaussian function Pf (a) = at a £ [0, 1/A^] 
and Pf (a) = 1 at a G [l/A^, . 

The asymptotic behavior at large a depends on the position of the pole closest to the real axis of the function 
in the s space. In the region of parameters (|B3p it is defined by the complex pole and has an oscillating character 
modulated by the exponential decay. The similar relation between " confining" (better to say screening) properties of 
propagator and complex singularities of the propagators in momentum space has been obtained in 71] for QED in 
three-dimensions for the fcrmion propagator and in |53j | for QCD in the instanton model for the gluon field strength 
correlator. ^ 

For the set of parameters y < yo the large a asymptotic of P (a) is purely exponential 

P (a) ~ exp{~KM^a). (B8) 

The function Pf (a) has the same asymptotics because in the s space it has the same dominating singularity. Through 
the use of the a-representation we avoid complexities typical for the momentum integral calculations. 



23 




FIG. 15: The propagator functions P (solid line) and Pf (dashed line) in a-representation. The a-representation of the 
propagator function P for "confining" set of parameters (A/q = 350 MeV, A = 815 MeV) is given by dotted line. 



APPENDIX C: QCD EVOLUTION 



Here we collect the formulas needed for the leading-order QCD evolution. The Gegenbauer coefficients evolve as 



aUQ')=<(Qo) 



^(g2)x(7;-7o)/(2/3o) 

a{Ql 



(CI) 
(C2) 



with 



r.3/2 



2n-|-3 



3(n + l)(n + 2) 

The anomalous dimensions are as follows (see for example in [t^): 



(C3) 
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^ AV 
In In 
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,0 I J. 



fe=i 
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